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c3 1. Introduction 



ABSTRACT 



Context. The dust produced by the first generation of stars will be a foreground to cosmic microwave background. 
Aims. In order to evaluate the effect of this early dust, we calculate the power spectrum of the dust emission anisotropies and compare 
it with the sensitivity limit of the Planck satellite. 

Methods. The spatial distribution of the dust is estimated through the distribution of dark matter. 

Results. At small angular scales (I > 1000) the dust signal is found to be noticeable with the Planck detector for certain values 
of dust lifetime and production rates. The dust signal is also compared to sensitivities of other instruments. The early dust emission 
anisotropies are finally compared to those of local dust and they are found to be similar in magnitude at mm wavelengths. 

Key words. Dust - CMB - Reionization - Power spectrum 
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The importance of the cosmic microwave background (CMB) 
as a cosmological tool has been demonstrated thoroughly dur- 
ing the last few years. It has been used to evaluate the age, the 
Hubble parameter, the baryon content, the flatness and the opti - 
cal depth of the reionization of the universe (iBennett et al.l2 003). 
It has also been used to set upper limits on the n on-Gaussianity 
of the primary fluctuations dKomatsu et al.ll2.Q03h, the Sun yaev- 
Zeldovich fluctuations from the first stars (lOh et aTl 2003), the 
primordial magnetic fields (Subramanian et al. 2003), the spa- 
tial curvature of the universe (|Efstathioul2003l) . the formation of 
population III stars (Cen 20 Q3|), and the neutrino masses using 
only WMAP data (Ichikawa et al. 2005) as well as combining 
them with other data (lHannes tad 2003). However, in order to in- 
terpret the CMB signal correctly, its foregrounds must also be 
well known. 

In this paper we focus on one particular aspect of the fore- 
ground of the CMB: the dust from the first generation of stars. It 
is here assumed that dust was created during the reionization pe- 
riod in the first generation of stars and was then ejected into the 
intergalactic medium (IGM). The dust is heated by the ionizing 
photons to a temperature slightly above Tqmb- The net effect on 
the CMB is a small monopole distortion of the CMB with a char- 
acteristic electromagnetic spectrum close to the CMB primary 
anisotropy (AT) sp ectrum times the frequen cy squared. This ef- 
fect was studied in lElfgren & Desertl d2004l) . 

Moreover, the dust also has a characteristic spatial dis- 
tribution, which could be used to identify its signal. The 
distribution gives rise to anisotropies in the dust emission, 
which can be measured with several current and future exper- 
iments. The objective of this paper is to determine this spa- 
tial distribution and its resulting anisotropies. Of particular in- 
terest is the Planck satellite mission, but other instruments 
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are a lso useful, like ALMA (IWoottenll2003h. BLAST dPevlinl 
120011) . BO LOCAM (LMT) a nd (CSO) dMauskopf et aljB oOO). 
MAMBO dGreve et alj|2004l) . SCUBA dBorvs et al.lll999l) . and 
SCUBA2 (tAudlev et al]l2004l) . 

The spatial distribution of the dust is estimated with the help 
of GallCS (Galaxies In Cosm ological Simulation s) A^-body sim- 
ulations of dark matter (DM) (Hat ton et al.l l2003). which are de- 
scribed in more detail in Sect. [2] The dust distribution is then 
combined with the intensity of the dust emission, and this is 
integrated along the line of sight. The resulting angular power 
spectrum is then computed as C{ and compared with the detec- 
tion limits of Planck. 

In the following, we assume a ACDM universe with D, to , = 
Q m + = 1, where Q,„ = Q/, + Sl DM = 0.133//i 2 , Q.b = 
0.0226/fr 2 , h = 0.7 2 and r e = 0.12, as advocated by WMAP 
dSpergel et al.l [2003). using WMAP data in combination with 
other CMB datasets and large-scale structure observations (2dF- 
GRS + Lyman a). 



2. Dark matter simulation 

The distribution of DM in the universe is calculated with the 
GallCS program. The cosmological N-body simulation we refer 
to throug hout this pa per is done with the parallel tree-code devel- 
oped by Ninin ( 1999). The power spectrum is set in agreement 
with lEke etalld 19961) : cr 8 = 0.88, and the DM-density field was 
calculated from z - 35.59 to z = 0, outputting 100 snapshots 
spaced logarithmically in the expansion factor. 

The basic principle of the simulation is to randomly dis- 
tribute a number of DM-particles A^ 3 with mass Mum in a box 
of size L 3 . Then, as time passes, the particles interact gravita- 
tionally, clumping together and forming structures. When there 
are at least 5 particles together, we call it a DM-clump. There are 
supposed to be no forces present other than gravitation, and the 
boundary conditions are assumed to be periodic. 
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In the simulation we set the side of the box of the simulation 
to L — lOOhr 1 Mpc and the number of particles to 256 3 , which 
implies a particle mass of ~ 5.51 x 10 9 hr x M Q . Furthermore, 
for the simulation of DM, the cosmological parameters were 
Q A = 2/3, Q„, = 1/3 and h = 2/3. The simulation of the DM 
was done before the results from WMAP were published, which 
explains the difference between these parameters and the val- 
ues used elsewhere in this paper, as stated in the introduction. 
Fortunately, the temporal distribution of the dust is independent 
of the value of h, which means that the impact of this small dis- 
crepancy is not important. Between the assumed initial dust for- 
mation at z ~ 15 and the end of this epoch in the universe at 
Z ~ 5, there are 51 snapshots. In each snapshot a friend-of-friend 
algorithm was used to identify virialized groups of at least 5 
DM-particles. For high resolutions, it is clear that the mass reso- 
lution is insufficient. Fortunately, the first 5-particle DM-clump 
appears at z = 14.7, while the bulk of the luminosity contribu- 
tion comes from z ^ 12.5. At z = 12.6 there are 19 clumps and at 
Z = 12.2 there are 45 clumps (and rapidly increasing in number) 
so it should be statistically sufficient. 

In order to make a correct large-scale prediction of the dis- 
tribution of the DM and therefore the dust, the size of the box 
would have to be of Hubble size, i.e., ~ 3000/!~' Mpc. However, 
for a given simulation time, increasing the size of the box and 
maintaining the same number of particles would mean that we 
lose out in mass resolution, which is not acceptable if we want 
to reproduce a fairly realistic scenario of the evolution of the 
universe. 

There is another way to achieve the desired size of the sim- 
ulation without losing out in detail or making huge simulations. 
This method is call ed MoMaF (Mock Map Facility) and is de- 
scribed in detail in Blaizot et al. (2005). The basic principle is 
to use the same box, but at different stages in time and thus a 
cone of the line of sight can be established. In order to avoid 
replication effects, the periodic box is randomly rotated for each 
time-step. This means that there will be a loss of correlation in- 
formation on the edges of the box, since these parts will be grav- 
itationally disconnected from the adjacent box. Fo rtunately, this 
loss w ill only be of the order of 10%, as shown in lBlaizot et alj 
(2005). For scales larger than the size of the box, there is obvi- 
ously no information whatsoever on correlation from the simu- 
lation. 

2.1. Validity of the simulation 

GallCS is a hybrid model for hierarchical galaxy formation stud- 
ies, combining the outputs of large cosmological N-body simu- 
lations with simple, semi-analytic recipes to describe the fate 
of the baryons within DM-halos. The simulations produce a de- 
tailed merging tree for the DM-halos, including complete knowl- 
edge of the statistical properties arising from the gravitational 
forces. 

The distribution of galaxies resultin g from this GallCS simu- 
lation has been compared with the 2dS ( Colless et alj|200ll) and 
the Sloan Digital Sky Survey dSzapudi et al.ll2001l) and found 
to be realistic on t he angular scales of 3' < 9 < 30', see 
Blaiz ot et alj ([2006). The discrepancy in the spatial correlation 
function for other values of 9 can be explained by the limits of 
the numerical simulation. Obviously, any information on scales 
larger than the size of the box (~ 45') is not reliable. Fortunately, 
the dust correlations increase at smaller angles, while the CMB 
and many other signals decrease. This means that our lack of 
information on angular scales 9 > 45' {( < 250) will not be 
important, as can be seen in Fig. [4] The model has also proven 
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Fig. 1. Intensity contribution from the dust per time-step z inte- 
grated over all frequencies. This model assumes that the fraction 
of metals produced in stars that end up as dust is fd = 0.3. The 
mean dust lifetime is a largly unknown parameter, and therefore 
three different values are explored, Af = 0.1, 1 and 10 Gyr. 



to give reasonable results for Lyman break galaxies at z — 3 
dBlaizot et al.i r2004). It is als o possible to model active galactic 
nuclei using the same model (Cat taneo et al.ll2005l) . 

Since it is possible to reproduce reasonable correlations from 
semi-analytic modeling of galaxy formation within this simula- 
tion at z — - 3, we hereafter attempt to do so at higher z values, 
when the early dust is produced. 



3. Model 

Since very little is known about the actual distribution of the dust 
throughout the universe at this time, we simply assume that the 
dust distribution follows the DM-distribution. We propose and 
explore two different ways for the dust distribution to follow the 
DM-distribution. The first is to let the dust be proportional to the 
DM-clumps, the second is to make a hydrodynamical smoothing 
of the DM-density field and set the dust density proportional to 
this density. In both cases we assume that 



Pdust(r,z) °c pDM(r,z), 



(1) 



where pbm represents either the clump method or the smooth- 
ing method density. The stars that produce the dust are proba- 
bly formed close to gravitational hot-spots in the universe and 
these spots are shaped by the DM. This means that the produc- 
tion of dust ca n be well approximate d by the DM distribution. 
According to lVenkatesan et all (120061) . the initial velocity of the 
dust is in the order of 10 5 m/s and the DM-simulations give halos 
with an escape velocity that is v e > 10 5 m/s, given that their mass 
is > 4.1 x 10 10 M . This means that most of the dust will stay 
near the halo and we will therefore focus on the clump method. 
The hydrodynamical smoothing case is included for reference 
only. 

In order to estimate the measured intensity, we need to cal- 
culate this distribution in terms of the intensity from the dust 
emission. The early dust is optically thin, a nd its intensity as 
a func tion of redshift has been calculated in lElfgren & Desert! 
(2004) and is shown in Fig.Q] This model assumes that the frac- 
tion of metals produced in stars that end up as dust is fd = 0.3. 
The mean dust lifetime is a largly unknown parameter and there- 
fore three different values are explored, Af = 0.1, 1, 10 Gyr. 
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In our present model, we put the spatial distribution of the 
dust intensity to 
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where dl(z)/dz is the dust intensity at redshift z as measured at 
z — and (pdm)(z) is the mean DM-density at redshift z. The 
MoMaF method (see Sect. [2J is then used to project the emitted 
intensity from the dust on a 45' x 45' patch along the line of 
sight. The contribution from each simulated box is added, and 
the integrated dust intensity is calculated. 

For z > 2.3, the time-steps are smaller than the size of the 
box and each box overlaps with the next box along the line of 
sight. However, for z < 2.3 the time-steps were simulated too 
far apart and when we pile the boxes, there will be a small part 
of the line of sight that will not be covered. Fortunately, this is 
of little consequence since the dust intensity is low at this time, 
and the gap is small. Each box (45' x 45') is divided into a grid 
according to the resolution that we wish to test. For Planck this 
means a grid that is 18x18 pixels 2 , for SCUBA 45x45 pixels 2 . 

To check the normalization of the resulting intensity image, 
we have calculated its £ dI x ,ylN 2 pix , where dl x>y is the observed 

intensity on pixel (x, y) and N 2 pjx is the number of pixels 2 , and 

found it to be equal to J dl(z)dz to within a few per cent. 



4. Results and discussion 

As described above, the MoMaF technique produces an image of 
the line of sight. This image represents the patch of the sky cov- 
ered by the box, 150 co-moving Mpc 2 , which translates to ~ 45 
arcmin 2 at z = 14.7. In order to avoid artifacts at the edges, the 
image is apodized, whereafter it is Fourier transformed into fre- 
quency space Pk- In order to convert this spectrum into spherical- 
harmonics correlation functions we apply the following transfor- 
mation: 



I = k2n/6, 



(3) 
(4) 



where 8 is the size in radians of the analyzed box. These Q are 
then calculated at a frequency v = 353 GHz, which is one of the 
nine P lanck frequency channels. As found in Elfgre n& Desert! 
(2004), the intensity is proportional to the frequency squared, 
which means that the power spectrum from the dust at a fre- 
quency v is 



Q(v) = Q(353 GHz).(^) 4 , 



(5) 



where Ct is given in terms of ^q^q ■ I n order to estimate an 
average power spectrum, 400 such images were generated and 
the C( were averaged over these. For comparison, we also tried 
to paste all these images together and calculate the Q for this 
(180x180 pixels 2 ) image. The result was similar to the average 
Ce . To validate our results, we also calculated the variance of the 
images and compared them with Ytt ^f^Qj and found them to 
be compatible. The resulting power spectra can be seen in Fig. [2] 
The lifetime of these dust particles is a largely unknown factor, 
and we plot three different lifetimes, 0.1, 1 and 10 Gyrs (fo r 
a more detailed discussion of dust lifetimes, see lDraind [T990). 
Furthermore, the dust intensity is proportional to the fraction of 
the formed metals that actually end up as dust, which we have 
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Fig. 2. Dust power spectrum in CMB thermodynamic units at 
353 GHz for a map 45' x 45' and Planck resolution 5' for three 
different lifetimes for the dust particles, 0.1, 1, and 10 Gyrs, with 
a solid, dashed and dotted line, respectively. The DM smoothing 
method for a dust lifetime of 1 Gyr is the dot-dot-dash line. 

assumed to be fj = 0.3. This means that the dust power spectrum 
is 



Ctifd) = C e (fd = 0.3) • (/ rf /0.3) 2 



(6) 



We note that there is only a small difference between dust life- 
times of 10 Gyrs and 1 Gyr, while the one at 0.1 Gyr is lower 
by a factor of four. The lowest curve in the figure represents the 
hydrodynamical smoothing method of distributing the dust for a 
dust lifetime of 1 Gyr. Naturally, it is lower than the correspond- 
ing Cc for the clump method, since the DM-clumps are much 
more grainy (especially early in history) than the smoothed DM 
field. The power spectra of the two methods differ by a factor of 
~ 10 but they do not have exactly the same form. 

The dust frequency spectrum is distinctly different from that 
of other sources in the same frequency range. In Fig. [3] we 
compare this spectrum with that of the primo ridal CMB AT 
aniso tropics and that of galactic dust, T — 17 K (Boulanger et al.l 
1996). In order to focus on the forms of the spectra, we normal- 
ize the three curves to one at v = 353 GHz. In case of a weak 
early dust signal, this frequency signature could help us identify 
the signal by component separation spectral methods. 

4.1. Detectability with Planck 

The Planck satellite^, due for launch in 2007, will have an angu- 
lar resolution of 30' -5' and will cover the whole sky. The Planck 
high frequency instrument (HFI) will measure the submillimeter 
sky at v = 100, 143, 217, 353, 545, and 857 GHz. We have cho- 
sen v = 353 GHz as our reference frequency. At higher frequen- 
cies, the galactic dust will become more of a problem, and at 
lower frequencies the CMB primary anisotropies will dominate. 

In order to test the detectability of the dust with Planck, we 
evaluate the sensitivity, 



v = o- c f m + aT tr 



X (E C MB + 



(7) 



" ^(2£+\)f cut L 

(iTegmarkl 1 19971) where f c u t — 0.8 is the fraction of the 
sky used, L is the bin-size, E C mb = €(€ + l)Cf MB /2n 



Homepage: http://www.rssd.esa.int/index.php ?project=Planck 
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Fig. 3. The form of the early dust spectrum compared to the form 
of galactic dust (with a temperature of 17 K) and the CMB along 
with indicators of the Planck HFI frequencies. The curves have 
been normalized to unity at 353 GHz. We see that the early dust 
has a special spectral signature. 
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Fig. 4. Comparison between dust power spectrum, the Planck er- 
ror limits and local dust (T = 17 K) at 353 GHz with binning 
500. The error limits (total noise) consist of two parts; the CMB 
cosmic variance, which dominates for small t, and the instru- 
ment noise, which dominates for high I. 



Tabl e 1. Parameters of the PL ANCK HFI detector properties 
dThe Planck collaborationll2005l) 



Frequency [GHz] 100 143 217 353 545 857 
FWHM['] 9.5 7.1 5.0 5.0 5.0 5.0 
sx L"Ks' /2 ] 32.0 21.0 32.3 99.0 990 45125 



(lLambda web-site: h ttp://lambda.g sfc.nasa.govl ll March 2 005) 
is the cosmic variance. The instrument error is 

p _ r - m „ A*+V _ f 4 *4 A-^ m 

tinstr - L ( - Jskx~ e * ■ , (8) 

In ' t bs In 

where / SJ t v = 1 is the fraction of the sky covered, Sx is the noise 
[//Ks 1/2 ], t obs = 14 ■ 30 • 24 • 3600 s is the observation time (14 
months), and cr b = FWHM/235 (FWHM is the full width at 
half maximum of the beam in radians). For Planck, the values of 
these parameters are given in Table Q] 

The resulting errors for a binning of L — 500 along with 
the dust power spectrum is plotted in Figs. |4]|5] In Fig. |H the 
frequency v = 353 GHz is fixed, while t is varied. We note that 
i ~ 1000 seems to be a good value to search for dust. At low I, 
the error due to the cosmic variance dominates, whereas at high 
t the instrument noise dominates. 

In Fig. [51 the i multipole binning center is fixed, and we show 
the electromagnetic spectrum of the primordial anisotropies and 
the early dust emission. The fourth point from the left in the 
figures corresponds to v = 353 GHz and gives the best signal 
over noise ratio. At low i the cosmic variance is important, at 
high €, the instrument noise dominates. 

Early dust may therefore produce a measurable disturbance 
in the primordial anisotropy angular power spectrum at high 
multipoles (in the Silk's damping wing). Although the primary 
contaminant to the CMB in the submillimeter domain is the in- 
terstellar dust emission, this new component vindicates the use 
of more than two frequencies to disentangle CMB anisotropies 
from submillimeter foregrounds. 

Component separation methods for the full range of Planck 
frequencies should be able to disentangle the CMB anisotropies 
from early dust, far-infrared background fluctuations, and galac- 
tic dust emission. The Sunyaev-Zel'dovich (SZ) effect also rises 



in the submillimeter range but the two central frequencies of the 
dust, 217 and 353 GHz, can be used for its identification since 
the SZ-effect drops more sharply at higher wave-lengths (the SZ- 
effect is practically null at 217 GHz, quite different from the 
early dust spec trum). The spectral sh ape of other foregrounds 
can be found in lAghanim et al.l (12000, Fig. 1) where we can also 
see that the region around I ~ 2000 is the most favorable for 
dust detection. 

4.2. Detectability with other instruments 

There are several other instr uments that mi ght be used t o de- 
tect t he early dust: ALMA (IWoottenl 120031 BLAST (iDevlinl 
120011) . BO LOCAM (LMT) a nd fCSOl jlMauskopf et alJl2000h . 
MAMBO dGreve et alj|2004 . SCUBA dBorvs et al.lll999h . and 
SCUBA2 (Aud lev et alJl2004h . 

Using Eqs. [7] and [8] we have estimated the sensibilities of 
these detectors. The result is presented in Table [2] Since all of 
these instruments operate on a small patch in the sky we use 
f sky = f cut = 10 ■ FOV, where FOV is the field of view of the 
instrument. The integration time was set to one hour and the 
noise per second, Sx, was calculated as sx = NEFD/ y/Nd e t, 
where Ndet is the number of detectors and NEFD is the noise 
equivalent flux density. The error was evaluated at the multi- 
pole moment t ~ 1 /FWHM (for BOLOCAM(LMT), I was 
set to 20,000), and we used a bin-size of L — I. Note that 

4-7TS fff+li 

ALMA is an interferometer and thus E jnstr = f S ky-r^ ' S 

'obs *£t 

Furthermore, the SCUBA2 array needs to be renormalized such 

that NEFD = NEFD x y/N de i/(FOV/(n/4 ■ {FWHM/60) 2 )). 
The resulting sensitivities cr'" s "' can be compared with the dust 
signal, as plotted in Fig. [6] As can be seen, BLAST, SCUBA 
and MAMBO are unable to detect the dust signal. However, 
BOLOCAM(LMT), ALMA, SCUBA2, and BOLOCAM(CSO) 
have good chances of detecting the radiation from the first dust. 
We also note that the curves are almost parabolic. In fact, for 
1, 000 < { < 100, 000 the curves can be fitted within ~ 10% as: 

tit + l)Cf l,s '/2n * 2.13 x 10~ 5 x £ Ln , At = 10 Gyr, 

{({ + l)C°'""/2n * 1.37 x 10~ 5 x £ 195 , At = 1 Gyr, (9) 

€(€ + l)Cf us, /2n ~ 4.02 x 10~ 5 x £ xm , At = 0.1 Gyr, 
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Fig. 5. Comparison between dust power spec- 
trum and Planck error limits at I = 1073, 
1358, 1731, and 2110 with binning 500. The 
error limits (total noise) consist of two parts; 
the CMB cosmic variance, which is constant 
at such ^K 2 -levels, and the instrument noise, 
which has the shape of an exponentialxf 2 . 



in units of fd^^E ■ The dependency in I, which is slightly dif- 
ferent from an uncorrelated noise ( 2 behavior, means that large- 
scale correlations cannot be neglected. They mix differently at 
different epochs, depending on the dust lifetime parameter. 



Table 2. Sensitivities, crf str , for different (current and future) de- 
tectors. NEFD = noise equivalent flux density, v is the operat- 
ing instrument frequency, Ndet — number of detectors, FWHM 
f Li ] ] width at half max, FOV = field of view in units of arcmin 2 , 
t = 1 1 FWHM, <t'" s "' is the instrument sensitivity in units of 
m K"CMB' The instrument sensitivity was calculated with Eqs.|7] 
&Husing t obs = lh, f sky = f cut = 10 • FOV and L-i. 
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5. Conclusions 

It seems that it is possible to detect the dust from the first gen- 
eration of stars with the Planck satellite on small angular scales 
(( > 1000). However, the detectability depends on the actual dis- 
tribution of dust in the early universe, and also, to a large extent, 
on the dust lifetime. 

The results are parametrized so that changing the frequency 
and the fraction of produced metals that become dust is only 
a matter of scaling the figures: Q oc (v/353 GHz) 4 and Q oc 
(fd/0.3) 2 . The spectrum of the early dust is compared to that 
of the primary CMB anisotropics, as well as that of the local 
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Fig. 6. Dust power spectrum as a function of multipole moment, 
{, for different dust lifetimes. Estimated detection limits for dif- 
ferent instruments are also included as small circles placed at 
6 ~ 1 1 FWHM. 



dust. The unique spectral signature of the early dust will help 
in disentangling it from the CMB and the different foregrounds 
(local dust and extragalactic far infrared background). 

The spatial signature of the early dust is found to have Q « 
constant as 10 - 10 -4 /"K^^g depending on the dust lifetime, 
At. Obviously, other signals that are correlated with the struc- 
tures will also show a similar behavior in the power spectrum. 
Notably, the near infrared background from primordial galaxies 
could be correlated with the early dust. 

The next generation of submillimeter instruments will be ad- 
equate to measure these early dust anisotropies at very small an- 
gular scales (€ > 2000). Our estimation shows that BOLOCAM, 
SCUBA2 and ALMA have a good prospect of finding the early 
dust. However, for these instruments, more detailed simulations 
are required in order to obtain a realistic DM and baryon distri- 
bution. A DM simulation on a smaller box, maybe L = 50/i _I 
for PLANCK and smaller still for ALMA, would improve the 
results on the relevant angular scales, £ > 1000. This also 
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means that the particles are smaller, giving a better level of de- 
tail. Furthermore, the distribution of dust relative to the DM can 
also be improved and it is even possible to include some semi- 
analytical results from the galaxy simulations in GallCS. 

Acknowledgements. Finally, we wish to express our gratitude towards Eric 
Hivon who wrote the code to analyze the simulated signal and thus produce 
the dust power spectrum. 
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